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If the observable universe is a braneworld of Randall-Sundrum 2 type, then particle interactions 
at high energies will produce 5-dimensional gravitons that escape into the bulk. As a result, the 
Weyl energy density on the brane does not behave like radiation in the early universe, but does 
so only later, in the low energy regime. Recently a simple model was proposed to describe this 
modification of the Randall-Sundrum 2 cosmology. We investigate the dynamics of this model, and 
find the exact solution of the field equations. We use a dynamical systems approach to analyze 
global features of the phase space of solutions. 
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I. INTRODUCTION 

String and M theory indicate that gravity may be a 
higher-dimensional interaction, while matter fields and 
gauge interactions are confined to a l-|-3-brane mov- 
ing in the higher-dimensional "bulk" spacetime. These 
ideas have led to relatively simple phenomenological 
models that allow one to explore some of the conse- 
quences of braneworld gravity. The second Randall- 
Sundrum (RS2) scenario is a five-dimensional Anti-de 
Sitter (AdSs) spacetime with an embedded Minkowski 
brane where matter fields are confined and Newtonian 
gravity is effectively reproduced at low energies. The RS2 
scenario is readily generalized to a Friedmann-Robertson- 
Walker (FRW) brane, and the bulk is Schwarzschild- 
AdSs. The modified Friedmann equation (for a flat uni- 
verse without cosmological constant) is 
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where A is the brane tension and C is the mass parameter 
of the bulk black hole. The tidal fleld of this black hole on 
the brane is felt as an effective radiative energy density, 
the so-called "dark radiation" . 

If the reheating and early radiation eras occur at suf- 
ficiently high energies, this simple picture is modified. 
High energy (p 3> A ^ (1 TeV)"*) particle interactions can 
produce 5D gravitons which are emitted into the bulk. 
This process will contribute to an increase in the dark 
radiation term, so that C will not be a constant in the 
early universe, but will grow until the energy scale drops 
below the threshold for 5D-graviton production. In the 
low-energy universe, C will have its asymptotic constant 
value, which is constrained by nucleosynthesis limits 0. 

In U , a simple model is constructed in which all the 5D 
gravitons are assumed to be emitted radially, and their 
backreaction on the bulk metric is assumed negligible, so 
that the bulk is a Vaidya-AdSs spacetime. (Earlier al- 
ternative models are developed in '3|.) Subsequently, the 
corrections to this model arising from non-radial emission 
have been investigated These corrections are very 
complicated, and we will confine ourselves to the simple 



radial model. We flnd the exact solution of the dynami- 
cal equations, and analyze the asymptotic behaviour and 
the global phase space of solutions. 



II. FIELD EQUATIONS 



The 5D field equations are 
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The projected field equations on the brane 0, 0] have 
three new terms from extra-dimensional gravity: 



G^u — ^^9iiv + + 6—5^^ — E^i, J'^v 1 (3) 

where k = + kI/\2 and ^ Ak|/6. The term 
Sfj,u is quadratic in T^j, and dominates at high energies 
(p > A); this term gives rise to the p^/A term in the 
modified Friedmann equation The five-dimensional 
Weyl tensor is felt on the brane via its projection, 5^^. 
This Weyl term is determined by the bulk metric, not 
by equations on the brane, and it gives rise to the C/a^ 
term in the modified Friedmann equation 

The term J-^j^u is also determined by the bulk metric, 
and arises from ^^^Tab, i-e., from 5D sources in the bulk 
other than the vacuum energy A5 , such as a bulk dilaton 
field. The presence of such sources modifies the energy- 
momentum conservation equations on the brane 



where is the unit normal to the brane. 
Radial radiation in the bulk is described by 
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where is a radial null vector and determines the 
energy density. A physical normalization of k^ is k^UA — 
1, where is the preferred 4- velocity, coinciding on the 
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brane with the 4-velocity of comoving observers. Then 
Eq. ||SJ) determines the form of J-^i, as 



(6) 



When the nuU radiation describes 5D gravitons produced 
in particle interactions in the radiation era (or in reheat- 
ing with radiative equation of state), then ^ 
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where a is a dimensionless constant. If ah degrees of 
freedom in the standard model are relativistic, then a ~ 
0.02. The brane energy-momentum tensor on the brane 
is assumed to be of the perfect-fluid type with respect to 



i.e. T„ 



For a FRW brane (which we assume to be spatially 
flat), the bulk metric has 4-dimensional plane symmetry, 
and then Eq. Q leads to the 5-dimensional Vaidya-AdSs 
metric 



(^)ds2 = -F{v,r)dv^ 
}'[v,r) — 



2drdv + r^dx^ 
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When C is independent of d, the metric reduces to the 
Schwarzschild-AdSs metric. The coordinate u is a null 
coordinate (corresponding to ingoing radial light rays). 
The brane trajectory is given by r = a(t) and v — v{t), 
where t is proper time for cosmic fluid particles, which 
move with 4-velocity u'*. The normalization of implies 



F{v,a) 
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For an expanding (outgoing) brane this is infinite at the 
apparent horizon {F{v, r) =0), so we assume that an ex- 
panding brane is always located in the region F{v,r) > 
(for Schwarzschild- AdSs this means outside the horizon) . 
See ^ for further discussion. 

The above equations lead to the modified Friedmann 
equation Q , and an equation governing the evolution of 
C (valid outside the apparent horizon) 



C ~ — aa p 



(11) 



The emission of 5D gravitons is at the expense of radia- 
tion energy density; by Eqs. Q-Q, 



p + AHp = — ^ ap^ . 



III. EXACT SOLUTION 
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We introduce the dimensionless variables 0: 

H = m, p 
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(Note that P = Q/n^X and = n'^d.) Then the high 
energy regime is given by p, > 1, and the dynamical 
equations Q, (fTT)l and ifT^ become 



a* 

^ = 2aa^p\l + p-H). 
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These equations imply the modified Raychaudhuri 
equation 



^ = -2i72-(a + 2)p^ 
dt 

which we use to decouple the system: 



-Ap-H) = 2[p-Hf 
dt 



This can be integrated to give 
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where to is an integration constant. Introducing this into 
Eq. IISI), we find 



dp 2p 



dt [t - to) 
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This equation is of Bernoulli type, with solution (restor- 
ing units) 



K^p 
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A{t - to){t - ti) ' 
where A is a dimensionless integration constant and 

,(4 + a) 
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Then, by Eq. 119|l the Hubble parameter can be written 
as 



H = 



1 



1 
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(23) 



a + 4:\t~ti 2(t-to) J ' 

and then the scale factor follows as 

a = ao{t - io)('+")/'^"+"^(t - ti)i/(4+") , (24) 

where oq is a constant (with dimension £^^^). Finally the 
expression for the dark radiation coefficient C is: 

^ at[A{A-8)it~t,)+M{A-2{4 + a)}] 

AA^{t~to)*/^'^+''Ht-ti)^/(^+'^') ' 

When ^ > (ti > to), the regions in which p > are: 
(i) t — ti < and t — to < i.e. — oo < t < to] and 
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(ii) t - ti > and t - to > 0, i.e. ti < t < oo. When 
A < (ti < tg), the p > region is: (iii) t — to < and 
t-h > 0, i.e. h <t < to. By Eq. ^ it follows that 
region (i) corresponds to contracting models that end in 
a big crunch att = tQ. The allowed part (where F{v, a) > 
0) of region (ii) corresponds to expanding models which 
start at t = > ii, where is the time at which the 
brane is just outside the apparent horizon (or at an initial 
singularity). The allowed part of region (iii) corresponds 
to models that expand from the apparent horizon t = , 
then collapse to a big crunch at i = tg. 

The ever-expanding models [region (ii)] are of most 
relevance. By Eq. (|24|l . it follows that near the singularity 
(t — > ii), which corresponds to the high-energy limit, 



ait) - (^-<i)i/(4+") 



(26) 



(we note that the solution is only valid near this limit if 
C is initially small) . This behaviour is different from that 
in the RS model, a(t) ~ f^/*, and the general relativistic 
one, a{t) ^ t^^^. Whereas p and H are always positive, 
the sign of C depends on the initial data, that is on the 
constant A. If < A < 8 then C is always negative. If 
8 < A < 8 -I- 2a, then C is initially negative and changes 
sign at i = <i + 4^[2(a -I- 4) - A]/yl(yl - 8), becoming 
positive forever. Finally, if A > 8 + 2a, then C is always 
positive. Normally one requires C > to avoid a naked 
singularity in the bulk. 

An interesting feature is the behaviour of the dark ra- 
diation U — C /a"* near the initial time . In general 
it has a minimum point at but grows quickly [this 
can be seen from Eqs. ((TH) and l(TT|l ]. For C > (i.e. 
A > 8 + 2a) the general case describes the growth of the 
black hole mass from a minimum initial state, when the 
brane is at the horizon. The solid line in Fig. ^illustrates 
this growth. 

However, for the special case A = 2(4 + 0;) the be- 
haviour is completely different: 
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This special case is the one identified by perturbative 
analysis in Ref. Q], but the general case was not given, 
and nor were the exact general solutions, Eqs. H21|) . 1)24(1 
and if^ . 

The special case corresponds toC = Oata = 0,i.e., 
the black hole is created from nothing. However, the cre- 
ation begins from a singularity (p(<i) ~ oo), so it is not 
clear what this means. The growth of C in this special 
case is illustrated by the dotted line in Fig. ^ 

In the low energy limit, the standard RS2 behaviour is 
recovered, as expected: 
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FIG. 1: Behaviour of the dark radiation coefficient C as a 
function oit~ti, in the general case ^4 > 8 + 2q: (solid curve) 
and the special case A = 8 -I- 2a (dashed curve). 



IV. DYNAMICAL SYSTEMS ANALYSIS 

Following 0, we use a dynamical systems approach 
to investigate the global properties of the phase space of 
solutions. Defining 



"^P ^ U C p 



the modified Friedmann equation H14|l takes the form 



ri[/ + = 1 . 



(30) 



In contrast to RS2 and general relativistic models, the 
phase space determined by the variables = {^Ip, ^u, ^) 
is not compact. This is because there are models in which 
ilu can change sign. We could restrict the phase space by 
the physical argument that > because C represents 
the mass of a black hole in the bulk. In that case, in the 
restricted phase space, the models in which C changes 
sign will be represented by trajectories that in the past 
cross the boundary of the restricted phase space. In all 
cases, Eq. (|3U|) implies that ^Ijj < 1. Whenever Ujj > 0, 
we must have < fip < 1 and —1 < uj < 1. 

It is convenient to rewrite the evolution equations with 
a new time parameter r, where 



dT ■' \H\ di ■ 



(31) 



Then 



H' ^ -s{l + q)H , (32) 
where s = sgn H and q is the deceleration parameter, 
1 (fa 



a2iJ2 dt-^ 



l + ia + 2)w' 



(33) 
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The evolution equations for are then 

n'p = sLu^lp [2{a + 2)uj - a] , (34) 

Lu' = suj [{a + 2)^2 ~-auj-2] . (36) 

Now the Friedmann equation l|3()|l is a constraint on 
the initial values of O. The equilibrium points of the 
dynamical system ft' ~ /(S^) determined by Eqs. H34|) - 
(|36|l are solutions of /(O*) = 0: 

= (0,0,1), (37) 

= (j/, l-iy,0). (39) 

Actually fig represents a set of equilibrium points 
parametrized by < < 1. The dynamical character 
of a given equilibrium point il* describes the dynamical 
behaviour in a small neighbourhood and is determined 
by the eigenvalues of the matrix ^fif ^* ■ 

The equilibrium point fij is a repeller for expanding 
models. Actually, since = 1 > 0, and considering only 
models with positive energy density, this point only exists 
for expanding models (s = 1 H > 0). From the fact 
that this equilibrium point is characterized by w = 1 it 
follows that a{t) oc and it represents the high 



energy regime [see Eq. (I2f)ll ]. which is different from the 
one found in RS2. Moreover, at this equilibrium point 
we have U = constant, so there is no flux of gravitons. 

The equilibrium point has ~ — 2/(a + 2) < 0, 
so that it only exists for contracting models (s ~ ^1)- 
At this point, a{t) cx {to - t)("+2)/[2(a+4)] and C{t) oc 
(to — i)~^/("+^). This equilibrium point, which is an at- 
tractor, is the analogous point to Jl^ for expanding mod- 
els. Therefore, there is an asymmetry between expand- 
ing and contracting models, contrary to what happens 
in RS2 and general relativistic models 0. This is es- 
sentially due to the fact that, unlike in RS2 and general 
relativity, the dynamical equations on the brane are not 
invariant under time reversal since there is energy loss to 
the bulk; specifically, the energy balance equation (|12|l is 
not invariant. This asymmetry can be seen in the planes 
Hp — u! and ri[/ — cj of the phase space shown in Figs. [21 
and 131 

The set of points rjg constitutes a line of equilibrium 
points, which are attractors (repellers) for expanding 
(contracting) models. At these points, a{t) oc i^/^, the 
same as in standard radiation models. 

The phase space is a non-compact three-dimensional 
space. For simplicity we only show the projections onto 
the three possible 2-planes: the flp — ui plane (Fig. [21), the 
ilfj — u! plane (Fig. [SJ, and the — flp plane (Fig.^. 
These planes describe completely the dynamics of the 
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FIG. 4: f2(7 — f2p plane of the phase space. The shaded region 
is where C > 0. 



models. In each of the figures, the (green) shaded region 
containing flow lines indicates where C is non-negative. 
In Figs, m and the (red) shaded region without flow 
lines is the region inside the horizon where our equations 
do not apply. Figures |21 and show the asymmetry be- 
tween expanding and contracting models, and how the 
dynamics is determined by the equilibrium points $7^ for 
expanding models (repeller) and for contracting mod- 
els (attractor). 

In Fig. ^ we can see the line of equilibrium points ff^ 
that extends to infinite values of Op and Qu- (Its projec- 
tion in the other two planes can be seen in Figs.[21and|2l. 
This plot describes the dynamics for contracting models 
only, where the point appears. The corresponding 



one for expanding models can be obtaining by replacing 
^2 by $7 J, that is, by moving the equilibrium point to 
the origin, and inverting the arrows. The three figures 
show the existence of trajectories on which flu, and hence 
C, changes sign, going from negative to positive values. 
Note that the contrary, to go from positive to negative 
values of C, is forbidden. Expanding models will always 
end up with a positive C. 



V. CONCLUSION 

We have investigated the dynamics of a RS2-type 
braneworld that is radiating 5D gravitons into the bulk 
during the high-energy radiation (or radiative reheating) 
era. The simplified model developed in 2] assumes that 
all gravitons are emitted radially, thus allowing the bulk 
spacetime to be described by the Vaidya-AdSs metric. 
We found the exact solution in this model to the field 
equations on the brane, as given in Eqs. H21 |l - H25|l . The 
Weyl radiation term C/a^ on the brane, due to the tidal 
effect of the bulk black hole, is only purely radiative (i.e., 
C — const) asymptotically, in the low-energy universe. 
At early times, C evolves due to graviton emission, with 
two cases that depend on the choice of an integration con- 
stant, as illustrated in Fig.^ Graviton emission leads to 
a loss of energy from the brane and a consequent asymme- 
try between expanding and collapsing models, that is not 
present in the non-radiating RS2 braneworld. This asym- 
metry, as well as the asymptotic behaviour of the models 
at early and late times, is illustrated in the phase planes 
that fully characterize the dynamics, shown in Figs.|2JSl 
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